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Abstract. We review recent applications of the real-space renormalization group 
(RG) approach to the integer quantum Hall (QH) transition. The RG approach, 
applied to the Chalker-Coddington network model, reproduces the critical distribu- 
tion of the power transmission coefficients, i.e., two-terminal conductances, Pc{G), 
with very high accuracy. The RG flow of P{G) at energies away from the transition 
yields a value of the critical exponent, vq = 2.39 ± O.Of, that agrees with most 
accurate large-size lattice simulations. Analyzing the evolution of the distribution 
of phases of the transmission coefflcients upon a step of the RG transformation, we 
obtain information about the energy- level statistics (ELS). From the fixed point of 
the RG transformation we extract a critical ELS. Away from the transition the ELS 
crosses over towards a Poisson distribution. Studying the scaling behavior of the 
ELS around the QH transition, we extract the critical exponent z^els = 2.37 ±0.02. 

The integer quantmn Hall (QH) transition is described well in terms of a 
delocalization-localization transition of the electronic wave functions. In con- 
trast to a usual metal-insulator transition (MIT), the QH transition is char- 
acterized by a single extended state located exactly at the center e = of each 
Landau band [1]. When approaching e = 0, the localization length ^ of the 
electron wavefunction diverges according to a power law e""^, where e defines 
the distance to the MIT for a suitable control parameter, e.g., the electron 
energy. On the theoretical side, the value of v has been extracted from vari- 
ous numerical simulations, e.g., v = 2.5 ±0.5 [2], 2.4 ±0.2 [3], 2.35 ±0.03 [4], 
and 2.39 ± 0.01 [5]. In experiments v k, 2.3 has been obtained, e.g., from the 
frequency [6] or the sample size [7] dependence of the critical behavior of the 
resistance in the transition region at strong magnetic field. 

Recently, a semianalytical description of the integer QH transition, based 
on the extension of the scaling ideas for the classical percolation [8] to the 
Chalker-Coddington (CC) model of the quantum percolation [2], has been de- 
veloped [9,10]. The key idea of this description, a real-space-renormalization 
group (RG) approach, is the following. Each RG step corresponds to a dou- 
bling of the system size. The RG transformation relates the conductance 
distribution of the sample at the next step to the conductance distribution 
at the previous step. The fixed point (FP) of this transformation, yields the 
distribution of the conductance, Pc{G), of a macroscopic sample at the QH 
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transition. This universal distribution describes the mesoscopic properties 

of a fully coherent QH sample. Analogously to the classical percolation [8], 
the correlation length exponent, u, was extracted from the RG procedure [5] 
using the fact that a slight shift of the initial distribution with respect to 
the FP distribution Pc{G) drives the system to the insulator upon rcnormal- 
ization. Then the rate of the shift of the distribution maximum determines 
the value of v. Remarkably, both Pc{G) and the critical exponent obtained 
within the RG approach [5,11,12] agree very well with the "exact" results of 
the large-scale simulations [3,4,13-15]. 

The goal of the present paper is threefold. First, we briefly review the basic 
ingredients that constitute the real-space RG method in the QH situation [5] . 
Second, we extend the RG approach to include the level statistics at the QH 
transition and apply a method analogous to the finite-size-corrections analysis 
to extract v from the energy-level statistics (ELS) obtained within the RG 
approach. This method yields v = 2.37 ± 0.02, which is even closer to the 
most precise large-scale simulations result v = 2.35 ± 0.03 [4] than the value 
V = 2.39±0.01 inferred from the conductance distribution [5]. This agreement 
is by no means trivial. Indeed, the original RG transformation [5] related 
the conductances, i.e., the absolute values of the transmission coefficients of 
the original and the doubled samples, while the phases of the transmission 
coefficients were assumed random and uncorrelated. In contrast, the level 
statistics at the transition corresponds to the FP in the distribution of these 
phases. Therefore, the success of the RG approach for conductances does not 
guarantee that it will be equally accurate quanUtatively for the level statistics. 
Third, we show that the RG structure employed in the present approach, 
which is constructed from 5 saddle points (SP), represents in many aspects 
the minimal model of the QH transition. A further reduction in the number 
of SP leads to less reliable results. 



1 The RG approach to the CC model 

Our RG approach to the QH transition [5,9,10] is based on the RG imit 
shown in Fig. 1. The unit is a fragment of the CC network consisting of 
five nodes. Each node, i, is characterized by the transmission coefficient ti, 
which is an amplitude to deflect an incoming electron along the link to the 
left. Analogously, the reflection coefficient = (1 — if )^^^ is the amplitude 
to deflect the incoming electron to the right. Doubling of the sample size 
corresponds to the replacement of the RG unit by a single node. The RG 
transformation expresses the transmission coefficient of this effective node, 
t' , through the transmission coefficients of the five constituting nodes [9] 

tih (rarsrie'*^ - l) + t2Ue^(^^+^^^ (rirarse-'*^ - l) 
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Fig. 1. Chalker-Coddington network on a squaxe lattice consisting of nodes (circles) 
and links (arrows). The RG unit used to construct the matrix (4) combines five 
nodes (full circles) by neglecting some connectivity (dashed circles). <Pi, . . . ,^4 are 
the phases acquired by an electron along the loops as indicated by the arrows. 
^1, . . . ,^4 represent wave function amplitudes, and the thin dashed lines illustrate 
the boundary conditions used for the computation of level statistics 

Here <Pj are the phases accumulated along the closed loops (Fig. 1). Within 
the RG approach to the conductance distribution, information about electron 
energy is incorporated only into the values of ti [5]. The energy dependence 
of phases, <Pj, is irrelevant; they are assumed completely random. Due to 
this randomness, the transmission coefficients, U, for a given energy, are also 
randomly distributed with a distribution function P{t). Then the transforma- 
tion (1) allows, upon averaging over ^j, to generate the next-step distribution 
P{t'). Therefore, within the RG scheme, a delocalized state corresponds to 
the FP distribution Pc (t) of the RG transformation. Due to the symmetry of 
the RG unit, it is obvious that the critical distribution, Pc{t^), of the power 
transmission coefficient, = G, which has the meaning of the two-terminal 
conductance, is symmetric with respect to ~ ^ as shown in Fig. 2. In other 
words, the RG transformation respects the duality between transmission and 
reflection. The critical distribution PdG) found in Refs. [5] and [9] agrees 
very well with the results of direct large-scale simulations. 

2 RG Approach to the ELS 

It has been realized long ago that, alongside with the change in the behavior 
of the eigenfunctions, a localization-delocalization transition manifests itself 
in the statistics of the energy levels. In particular, as the energy is swept 
across the mobility edge, the shape of the ELS crosses over from the Wigner- 
Dyson distribution, corresponding to the appropriate universality class, to 
the Poisson distribution. Moreover, finite-size corrections to the critical ELS 
close to the mobility edge allow to determine the value of the correlation 
length exponent [16], thus avoiding an actual analysis of the spatial extent 
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Fig. 2. The critical distribution of the conductance -Pc(G) at the QH transition 
obtained using the 5 SP RG unit (dashed line). The dotted line denotes a 4 SP 
RG unit as discussed in Sect. 4. The latter distribution clearly deviates from the 
expected symmetry with respect to G = 0.5 



of the wave functions. For this reason, the ELS constitutes an alternative 
to the MacKinnon-Kramer [17 20] and to the transmission-matrix [21,22] 
approaches to the numerical study of localization. 

2.1 Derivation of the Network Operator for the RG 

As has been shown by Fertig [23], energy levels of an 2D CC network can 
be computed from the energy dependence of the so called network operator 
U{E). U is constructed similar to the system of equations for obtaining the 
transmission coefRcient t' of the RG unit as presented in Eq. (1). Every SP of 
the network contributes two scattering equations. Each of them describes the 
amplitude of one outgoing channel using the amplitudes of the two incoming 
channels weighted by the transmission and reflection coeflicients t and r in 
which also the random phase <P of the links between SP's can be incorporated. 
When comparing to the calculation of the transmission coefficient t' an es- 
sential difference has to be taken into account. Energy levels are defined only 
in a closed system which requires to apply appropriate boundary conditions. 
The energy dependence of U (E) enters trough the energy dependence of the 
ti{E) of the SP's, as well as the energy dependence of the phases ^j{E) of 
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the links. Considering the vector ^ of wave amphtudes on the Hnks of the 
network, the eigenenergies can now be obtained from the stationary condition 

U{Ep = ^. (2) 

Nontrivial solutions exist only for discrete energies Ek, which coincide with 
the eigenenergies of the system [23]. The evaluation of the E^'s according to 
Eq. (2) is numerically very expensive. For that reason a simplification was 
proposed [24]. Instead of solving the real eigenvalue problem, calculating a 
spectrum of quasienergies co is suggested following from 

U{E)^i = e^'^^'^^i. (3) 

For fixed energy E the uji arc expected to obey the same statistics as the 
real eigenenergies [24] . This approach makes is perfectly suited for large-size 
numerical simulations, e.g. studying 50 x 50 SP networks. 

In order to combine the above algorithm with the RG iteration, in which 
a rather small unit of SP's is considered, we first "close" the RG unit at 
each RG step in order to discretize the energy levels as shown in Fig. 1 
with dashed lines. For a given closed RG unit with a fixed set of t^- values 
at the nodes, the positions of the energy levels are determined by the energy 
dependences, $j{E), of the four phases along the loops. These phases change 
by ^ TT within a very narrow energy interval, inversely proportional to the 
sample size. Within this interval the change of the transmission coefficients 
is negligibly small. A closed RG unit in Fig. 1 contains 10 links, and, thus, 
it is described by 10 amplitudes. Each link is characterized by an individual 
phase. On the other hand, it is obvious that the energy levels are determined 
only by the phases along the loops. One way to derive U is to combine the 
individual phases into phases <Pj connected to the four inner loops of the 
unit and to exclude from the original system of 10 equations all amplitudes 
except the "boundary" amplitudes (Fig. 1). The network operator for the 
remaining four amplitudes is a 4 x 4 matrix {Unm} with elements 
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(4) 



which can be substituted in Eq. (3). Then the energy levels, Ek, of the closed 
RG unit including phases '^j{E) = $j{Ek); are the energies for which one of 
the four eigenvalues of the matrix U is equal to one. Thus, the calculation of 
the energy levels reduces to a diagonalization of the 4x4 matrix. 
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The crucial step now is the choice of the energy dependence $j{E). If 
each loop in Fig. 1 is viewed as a closed equipotential as it is the case for the 
first step of the RG procedure [2] , then <l>j {E) is a true magnetic phase and 
changes linearly with energy with a slope governed by the actual potential 
profile, which, in turn, determines the drift velocity. Thus 

^o{E)=^o,j + 2n^, (5) 

where a random part, ^?o,j, is uniformly distributed within [0, 2it\, and 27r/sj 
is a random slope. Here the coefficient Sj acts as an initial level spacing 
connected to the loop j of the RG unit by defining a periodicity of the cor- 
responding phase. Strictly speaking, the dependence (5) applies only for the 
first RG step. At each following step, n > 1, ^j{E) is a complicated function 
of E which carries information about all energy scales at previous steps. How- 
ever, in the spirit of the RG approach, one can assume that ^j{E) can still be 
linearized within a relevant energy interval. The conventional RG approach 
suggests that different scales in real space can be decoupled. Linearization of 
Eq. (5) implies a similar decoupling in energy space. 

With ^j{E) given by Eq. (5), the statistics of energy levels determined 
by the matrix equation (3) is obtained by averaging over the random initial 
phases <l>oj and values ti chosen randomly according to a distribution P{t). 
For every realization the levels Ek are computed from the solutions Lo{Ek) = 
of Eq. (3) yielding 3 level spacings as illustrated in Fig. 3. Thus the situation 
is comparable with estimating the true random matrix ensemble distribution 
functions from small, say, 2x2 matrices only [25,26]. Within the RG approach, 
the slopes sj as in Eq. (5) determine the level spacings at the first step. 
They are randomly distributed with a distribution function Po{s). Subsequent 
averaging over many realizations yields the ELS, -Pi(s), at the second step. 
Then the key element of the RG procedure, as applied to the level statistics, 
is using Pi(s) as a distribution of slopes in Eq. (5). This leads to the next-step 
ELS and so on. 

The approach of this work relies on the real cigcnencrgies of the RG unit. 
The simpler computation of the spectrum of quasienergies adopted in large- 
scale simulations within the CC model [24,27] cannot be applied since the 
energy dependence of phases in the elements of the matrix is neglected and 
only the random contributions, ^oj, are kept. Nevertheless it is instructive 
to compare the two procedures as presented in Fig. 3. The figure shows the 
dependence of the 4 quasienergies LUk on the energy E calculated for two 
single sample RG units, with ti chosen from the critical distribution Pc{t). 
The energy dependence of the phases was chosen from the ELS of the 
unitary random matrix ensemble (GUE) according to Eq. (5). It is seen that 
the dependences (jj{E) range from remarkably linear and almost parallel (Fig. 
3a) to strongly nonlinear (Fig. 3b). 
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Fig. 3. Energy dependence of the quasieigenenergies ui for two sample configura- 
tions. Instead of using the quasispectrum obtained from u>i{E = 0) (O) we calculate 
the real eigenenergies according to u!{Eu) = (□). Different line styles distinguish 

different uJi{E). Wc emphasize that the observed behavior varies from sample to 
sample between remarkably linear (a) and strongly nonlinear (b) 



2.2 The Shape of the ELS at the QH Transition 

First, let us turn our attention to the shape of the ELS at the QH transition. 
As starting distribution Pq{s) of the RG iteration, wc choose the ELS of 
GUE, since previous simulations [24,28] indicate that the critical ELS is close 
to GUE. According to Po{s), Sj is drawn randomly and 'Pj, j = 1, . . . , 4 is set 
as in Eq. (5). For the transmission coefficients of the SP the FP distribution 
Pc{t), obtained in Sect. 1, is used as initial distribution Po{t). And from Po{t), 
the 5 ti, i = 1,...,5, are selected. As in Sect. 1 the RG transformation (1) 
is used to compute 10^ super-transmission coefficients t' . The accumulated 
distribution Pi(i') is again discretized in at least 1000 bins, such that the bin 
width is typically 0.001 for the interval t e [0, 1]. Pi{t') is then smoothed by 
a Savitzky-Golay filter [29] in order to decrease statistical fluctuations. By 
finding solutions u){Ek) = of Eq. (3) the new ELS Pi{s') is constructed 
from the "unfolded" energy-level spacings = {Em+i — Em)/ A, where 
m = 1, 2, 3, Sfe+i > Ek and the mean spacing A = [E^ — Ei)/3. Due to the 
"unfolding" [30] with A, the average spacing is set to one for each sample 
and in each RG-iteration step spacing data of 2 x 10^ super-SP's can be 
superimposed. The resulting ELS is discretized in bins with largest width 
0.01. In the following iteration step the procedure is repeated using the Pi's as 
initial distributions. Convergence of the iteration process is assumed when the 
mean-square deviations of both distributions Pn{i) and Pn{s) deviate by less 
than lO"** from predecessors P„_i(t) and P„_i(s). Once the (unstable) FP 
has been reached, the P„'s should in principle remain unchanged during all 
further RG iterations. Our simulations show [5] that unavoidable numerical 
inaccuracies sum up within several further iterations and lead to a drift away 
from the FP. In order to stabilize our calculation, we therefore use in every 
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RG step instead of P„(t) the FP distribution Pc{t). This trick does not alter 
the resuhs but speeds up the convergence of the RG for Pc{s) considerably. 

This now enables us to determine the critical ELS Pc{s). The RG iteration 
converges rather quickly after only 2 — 3 RG steps. The resulting Pc{s) is 
shown in Fig. 4 together with the ELS for GUE. Pc{s) exhibits the expected 
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Fig. 4. FP distributions -Pc(s) obtained from the spectrum of u)i{E = 0) and from 
the RG approach using the real eigenenergies Ek in comparison to the ELS for 
GUE. As in all other graphs P{s) is shown in units of the mean level spacing A 

features, namely, level repulsion for small s and a long tail at large s, but the 
overall shape of Pc{s) differs noticeably from GUE. In the previous large-size 
lattice simulations [24,28] the obtained critical ELS was much closer to GUE 
than Pc{s) in Fig. 4. This fact, however, does not reflect on the accuracy 
of the RG approach. Indeed, as it was demonstrated recently, the critical 
ELS although being system size independent nevertheless depends on 
the geometry of the samples [31] and on the specific choice of boundary 
conditions [32,33]. Sensitivity to the boundary conditions does not affect the 
asymptotics of the critical distribution, but rather manifests itself in the 
shape of the "body" of the ELS. Recall now that the boundary conditions 
which have been imposed to calculate the energy levels (dashed lines in Fig. 
1) are non-periodic. 

As mentioned above, the critical ELS has also been computed previously 
by diagonalizing U (0) and studying the distribution of quasienergies. In Fig. 
4 the result of this procedure using the present RG approach is shown. It 
appears that the resulting distribution is almost identical to Pcis). This ob- 




Fig. 5. Loft: FP Pc(s) for small s in agreement with the predicted behavior. 
Due to the log-log plot errors are shown in the upper direction only. Right: The 
large s tail of Pc{s) compared with fits according to the predictions of Ref. [16] 
(lines). The interval used for fitting is indicated by the bars close to the lower axis. 
For clarity errors are shown in upper direction and for s/A = 1.5, 2.0, 2.5, 3.0 only. 
For s/A < 2.4, only every 5th data point is drawn by a symbol 



servation is highly non-trivial, since, as follows from Fig. 3, there is no simple 
relation between the energies and quasienergies. 

2.3 Small and Large s Behavior 

As we have seen, the general shape of the critical ELS is not universal. How- 
ever, the small-s behavior of Pc{s) must be the same as for GUE, namely 
Pc{s) oc s^. This is because delocalization at the QH transition implies level 
repulsion [16,24,27,28,34-42]. In Fig. 5 we show that this is also true for the 
RG approach. The given error bars of our numerical data are standard de- 
viations computed from a statistical average of 100 FP distributions each 
obtained for different random sets of fj's and $j's within the RG unit. In 
general, within the RG approach, the ,s^-asymptotics of P{s) is most natural. 
This is because the levels are found from diagonalization of the 4x4 unitary 
matrix with absolute values of elements widely distributed between and 1. 
The right form of the large- .s tail of P{s) is Poissonian, Pc{s) oc exp(— &,s) [16] 
for s > 3 A. The data has a high accuracy only for s/zi < 2.5. For such s, the 
distribution Pc{s) does not yet reach its large-s tail and the fit parameters 
shown in Fig. 5 depend largely on the s-interval chosen. 

3 Scaling Results for the ELS 

3.1 Finite-Size Scaling at the QH Transition 

The critical exponent i' of the QH transition governs the divergence of the 
correlation length as a function of the control parameter zq, i.e. 

^oo{zo) OC \zo - Zcl'", (6) 
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where is the critical value. For the QH transition v k. 2.35 has been 

calciilatcd by a variety of numerical methods [2-4] and is in agreement with 
the experimental estimates u w 2.3 [6,43,44]. As presented in Sect. 1 the 
RG approach for the conductance distribution yields a rather accurate value 
V = 2.39 ± 0.01. In order to extract v from the ELS the one-parameter- 
scaling hypothesis [45] is employed. This approach describes the rescaling of 
a quantity a{N;{zi}) — depending on (external) system parameters {zi} 
and the system size N — onto a single curve by using a scaling function / 



Since Eq. (6), as indicated by "oo", holds only in the limit of infinite system 
size, we now use the scaling assumption to extrapolate f to N ^ oo from the 
finite-size results of the computations. The knowledge about / and ^cxa then 
allows to derive the value of u. 

We use the natural parametrization t = {e^ + [9], such that z can 

be identified with a dimensionless electron energy. The universal conductance 
distribution at the transition, Pc{G), corresponds to a distribution Qc{z) [5] 
which is symmetric with respect to ^; = and has a shape close to a Gaussian. 
The RG procedure for the conductance distribution converges and yields 
Qc{z) only if the initial distribution is an even function of z. This suggests to 
choose as a control parameter in Eq. (7), the position Zq of the maximum of 
the function Q{z). The meaning of zq is an electron energy measured from the 
center of the Landau band. The fact that the QH transition is infinitely sharp 
implies that for any zo 7^ 0, the RG procedure drives the initial distribution 
Q{z — zo) towards an insulator, either with complete transmission of the 
network nodes (for zq > 0) or with complete reflection of the nodes (for 
zo < 0). 

3.2 Scaling for ap and ai 

In principle, one is now free to choose for the finite-size scaling analysis (FSS) 
any characteristic quantity a{N; Zq) constructed from the ELS which has a 
systematic dependence on system size N for zq ^ while being constant 
at the transition zq = 0. Because of the large number of possible choices 
[16,28,46-49] a restriction to two appropriate quantities is made which are 
obtained by integration of the ELS and have already been successfully used 
in Refs. [46,50], namely 



with I{s) = P{s')d.s' . The integration limit is chosen as sq — 1-4 which 
approximates the common crossing point [46] of all ELS curves as can be seen 
in Fig. 6. Thus P(so) is independent of the distance jz — Zcj to the critical 
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and 
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Fig. 6. RG of the ELS used for the computation of v. The dotted lines corresponds 

to the first 9 RG iterations with an initial distribution Po shifted to the metallic 
regime [zq = 0.1) while the thin full lines represent results for a shift toward 
localization {zo = —0.1). Within the RG procedure the ELS moves away from the 
FP as indicated by the arrows. At s/A « 1.4 the curves cross at the same point - 
a feature we exploit when deriving a scaling quantity from the ELS 

point and the system size magnification TV. Since ai^p(A^, zq) is analytical for 
finite N, one can expand the scaling function / at the critical point. The first 
order approximation yields 

a{N,ZQ)r^a{N,z,)+a\zo-z,\N^''' (9) 

where a is a coefficient. For our calculation we use higher order expansions 
[51] expanding / twice, first, in terms of Chebyshev polynomials of order 

and, second, as Taylor expansion with terms |zo — -Zcl in the order Oz- 
This procedure allows to describe deviations from linearity in |zo — Zc| at the 
transition. Contributions from an irrelevant scaling variable can be neglected 
since the transition point zq = is known. In Fig. 7 the resulting fits for ap 
at the transition are shown. 

The fits are chosen in a way such that the total number of parameters 
is kept at a minimal value and the fit agrees well with the numerical data. 
The corresponding scaling curves for ap are displayed in Fig. 7. In the plots 
the two branches for complete reflection (zq < 0) and complete transmission 
{zq > 0) can be distinguished clearly. In order to estimate the error of the 
fitting procedure the results for v obtained by different orders O^, and of 
the expansion, system sizes N, and regions around the transition are com- 
pared. A part of the over 100 fit results together with the standard deviation 
of the fit are given in Ref. [52]. The value of u is calculated as average of all 
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Fig. 7. Left: Behavior of ap at the Qlf transition as results of the RG of the ELS. 
Data are shown for RG iterations n = 1, ... ,9 corresponding to effective system 
sizes A'^ = 2" = 2, . . . , 512. Full lines indicate the functional dependence according 
to FSS using the minimization with = 2 and Oz = 3. Right: FSS curves 
resulting from the fit of our data shown in Fig. 7. Different symbols correspond 
to different effective system sizes N = 2" . The data points collapse onto a single 
curve indicating the validity of the scaling approach 



individual fits where the resulting error of v was smaller than 0.02 resulting 
in 1/ = 2.37 ± 0.02. This is in excellent agreement with the previously quoted 
results [2-5]. 

4 Test of Different SP Unit 

Apparently, the quality of the RG approach crucially depends on the choice 
the RG unit. For the construction of a proper chosen RG unit two conflicting 
aspects have to be considered, (i) With the size of the RG unit also the 
accuracy of the RG approach increases since the RG unit can preserve more 
of the connectivity of the original network, (ii) As a consequence of larger 
RG units the computational eS'ort for solving the scattering problem rises, 
especially in the case where an analytic solution, as Eq. (1), is not attained. 
Because of these reasons building an RG imit is an optimization problem 
depending mainly on the computational resources available. As mentioned 
in the previous Section larger RG units were already studied in [11,12,53]. 
In these works the authors could not benefit from an analytic solution and 
achieve only a similar and less accurate statistics in comparison with the 
results presented here. In this Section the opposite case is studied using a 
small RG unit proposed in Rcf. [54] in the context of the Hall resistivity. 

The super-SP now consists only of 4 SP's as shown in Fig. 8. It resembles 
the 5 SP's unit (Fig. 1) used previously leaving out the SP in the middle of 
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Fig. 8. RG unit constructed from 4 SP's indicated by full circles. Some connectivity 
is neglected (dotted circles). The phases $j are accumulated by the electron motion 
(arrows) along contours of the energy potential 



the structure. Again the scattering problem can be formulated as a system 
of now 8 equations which is solved analytically 

(l-r2r3e**3)(i_rir4e**i)+*i*2i3t4e'*2 ' ^ ' 

The result can be verified using Eq. (1) after setting = and ra = 1, 
joining the phases <l>i and ^4 and renumbering the indices. 

The RG transformation (10) is then applied within the RG approach 
analogously to the 5 SP unit. First the FP distribution Pc{G) is obtained. A 
comparison of Pc{G) for both RG units is shown in Fig. 2. In the case of the 5 
SP unit the FP distribution Pc{G) exhibits a flat minimum around G = 0.5, 
and sharp peaks close to G = and G = 1. It is symmetric with respect to 
G « 0.5. The 4 SP unit yields differing results. While Pc(G) is stih rather 
flat it is clearly asymmetric, which already indicates that the 4 SP unit can 
not describe all of the underlying symmetry of the CC network. 

The Pc{G) for the 4 SP unit is then used in the calculation of the critical 
exponent v to construct the shifted initial distributions Qo{z)- The behavior 
of 1/ as function of n for the 4 and 5 SP RG units is demonstrates in Fig. 
9. Both curves approach convergence monotonously from larger values of u. 
During all iteration steps, u for the 4 SP differs from the 5 SP result by an 
almost constant positive shift. After 8 iterations, which equals an increase 
of system size by a factor of 256, one finds i^ssp = 2.39 ± 0.01 and V4sp = 
2.74 ±0.02. The error describes a confidence interval of 95% as obtained from 
the fit to a linear behavior. The result for i/4sp deviates clearly from the 
five SP result and also from the values obtained by other methods [2-4]. In 
addition to these findings also the discussion in Sect. 1 indicates that the 
4 SP RG unit fails to describe the critical properties at the QH transition 
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256 



Fig. 9. The critical exponent v as function of the effective system size iV = 2" 
for 4 SP (dotted hne) and 5 SP unit (dashed Une). Inset: Maximum Zraax. of Q{z) 
vs. initial shift for 8 RG iterations (symbols) using 4 SP. Dashed lines indicate 
linear fits 



correctly. This fact underlines again the importance of the RG unit for a 
successful application of the RG approach. 



5 Conclusions 

The version of the network model [55] that has been most widely studied 
in the context of the QH effect, is the CC model [2], describing the electron 
motion in a disordered system in a strong magnetic field limit. The fact that 
the RG approach, within which the correlations between different scales arc 
neglected, describes the results of the large-scale simulations of the CC model 
so accurately, indicates that only a few spatial correlations within each scale 
are responsible for the critical characteristics of the quantum Hall transition. 
More precisely, the structure of the eigenstates of a macroscopic sample at the 
transition can be predicted from the analysis of a single RG unit consisting of 
only five nodes. Further applications of this approach to the computation of 
the Hall resistance and the plateau-to-insulator transition shall be published 
elsewhere. 

We thank B. Huckestein, M.E. Raikh, M. Schreiber, and U. Ziilicke for 
stimulating discussions. This work was supported by the DFG within SFB393 
and the priority research program on quantum Hall systems. Further support 
was provided by a DAAD-NSF collaborative research grant INT-0003710. 
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